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Abstract 

CO 

We study the potential of a very long baseline neutrino oscillation experiment 
with a neutrino factory and a large segmented water-Cerenkov calorimeter detector in 
resolving the degeneracies in the neutrino oscillation parameters; the sign of the larger 
Oh. mass-squared difference 6m1 3 , the sign of \U^ 3 \ 2 (= sin 2 # ATM ) — 1/2, and a possible 

two-fold ambiguity in the determination of the CP phase <5 MNS . We find that the above 
problems can be resolved even if the particle charges are not measured. The following 
results are obtained in our exploratory study for a neutrino factory which delivers 
10 21 decaying /i + and fi~ at 10 GeV and a 100 kton detector which is placed 2,100 km 
away and is capable of measuring the event energy and distinguishing eP 1 from fi^, 
but not their charges. The sign of <5m 2 3 can be determined for 4|[/e 3 | 2 (l — | C^e3 1 2 ) = 
sin 2 26» RCT ^0.008. That of sin 2 <9 ATM - 1/2 can be resolved for sin 2 26» ATM = 0.96 
when sin 2 29 RCT J> 0.06. The CP-violating phase <5 MNS can be uniquely constrained 
for sin 2 2# RCT 0.02 if its true value is around 90° or 270°, while it can be constrained 



RCT 

for sin 2 2# RCT J> 0.03 if its true value is around 0° or 180° 
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Atmospheric neutrino observation at Super-Kamiokande (SK) pQ and the K2K exper- 
iment El established v., oscillation with a close-to-maximal mixing sin 2 29 > 0.9 and 

' — S At ° ATM 

5m 2 ATM ~ (1.3 - 3.0) x 10~ 3 eV 2 . The MSW large-mixing-angle (LMA) solution ij 
of the solar-neutrino deficit problem U\ has been established by KamLAND jH] and 
by the improved measurement by SNO jjj. The best fit values are sin 2 2# goL = 0.82 and 
Sm 2 QL = 7.1 x 10 -5 eV 2 . The CHOOZ 9] and Palo Verde ^U] reactor experiments give 
upper bounds on the third mixing angle of the three-neutrino model; sin 2 29 RCT ^0.1 (0.2) 
for 5m 2 ~ 3.0 (2.0) x 10~ 3 eV 2 . 

ATM v / 

In the three-neutrino model, the present neutrino oscillation experiments constrain U e 2, 
U e 3, and Ufxa matrix elements of the lepton-flavor mixing matrix U MNS (Maki-Nakagawa- 
Sakata (MNS) [TT]) : 



U, 3 \ 2 = sin 2 6 atm = (l±Jl- sin 2 26 ATM / 2 , (la 



U e2 \ 2 = [1 - \U e3 \ z - d (1 - \U e3 \ z )~ -sm 2 29 SOL \ / 2, (lb) 



U e3 \ 2 = 1- v /l-sin 2 2^ CT /2, (lc) 



and the mass-squared differences, Sm 2 QL = \5ml 2 \ <C \Sm1 3 \ = 5^ 2 TM , where Smjj = 
m 2 — m 2 . The matrix elements U e 2 and are taken to be real and non- negative while Z/ e3 
is a complex number in our convention [T2"l ll3j. The CP- violating phase <5 MNg = —arg{U e3 ), 
is unconstrained. Note that the solution for U e 2 follows from our convention JJ e \ > U e2 
[12], which defines the mass-eigenstate v v 

The present experiments allow several degeneracies in the neutrino oscillation parame- 
ters. First is the sign of the larger mass-squared difference, bm\ z = m 2 — m 2 = ±5m 2 ATM . 
The sign of the smaller mass-squared difference is determined as 8m\ 2 = m 2 — m 2 = 5m 2 ql 
by the MSW matter effects in the sun. The mass hierarchy 5m 2 3 = 5m 2 ATM is called 'normal' 
and Sm 2 3 = — Sm 2 ATM is called 'inverted'. The second degeneracy is in the mixing angle 
9 ATM . Although the present experiments constrain sin 2 2# ATM ~ 1, sin 2 9 has two-fold 



ambiguity, sin 2 ^ ATM — 1/2 = ±y 1 — sin 2 29 atm /2. Furthermore, there can appear two-fold 
ambiguity in the future determination of the <5 TO phase, between <S and 180° — S^ ma . 

° J MNS 1 ' MNS MNS 

Implications of the above degeneracy problems and their partial resolutions have been 
proposed in (2H|. 

The proposed long-baseline (LBL) neutrino-oscillation experiments [THITH] will achieve 
the precision measurements of Sm 2 ATM and sin 2 29 ATM by using conventional neutrino beams, 
which are made from decays of tt and K that are produced by high-energy proton beams. 
However, lifting the above eight-fold degeneracy in the neutrino oscillation parameters 
could remain as the target for the next generation neutrino experiments. 

In Ref. [16 , we studied the possibility of measuring 5 in the LBL experiment with 
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J-PARC (Japan Proton Accelerator Complex) at Tokai Village and a Megaton-level 
water- Cerenkov detector, Hyper-Kamiokande (HK) |18j . The distance between Tokai Vil- 
lage and HK is about 295 km. There we found that it is relatively easy to distinguish 
between 5 „ = 90° and 270° cases, but that it is difficult to distinguish between 5 „ = 0° 

MNS ° MNS 

and 180° cases. Those studies, however, assumed that the neutrino mass hierarchy is known 
to be normal (<5m?o = 5m 2 ) and also we assumed sin 2 9 = 0.5 as an input. 

V 16 ATM' ATM r 

In Fig^a), we show the allowed regions of Tokai-to-HK experiment when 5 MNS = 90° 
and the normal hierarchy is assumed. CP-violation can be established at better than 3a 
level if the neutrino mass hierarchy is known. However, the region encircled by thin lines 



are allowed if the inverted hierarchy (5m 2 3 



-5m 2 



is assumed in the analysis of the 



same data. CP conservation is preferred in the latter case. In Fig^b), we show the allowed 
region of the same experiment when sin 2 #* rue = 0.45, 0.35. Mirror solutions at around 



sin 9 



2 /imirror 



0.55, 0.65, respectively, are clearly seen* 
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Figure 1: The regions allowed by future LBL experiment between J-PARC at Tokai and 
HK in the plane of <5 MNS and sin 2 29 BCT (a), and sin 2 9 ATM and sin 2 29 RCT (b). The assumed 
experimental conditions are 0.8 Mton-year for OAB(2°), 2 Mton-year for OAB(3°), 
and 3.2 Mton-year for OAB(2°). x 2 fittings are performed by assuming the normal and 
the inverted hierarchies in (a) and the normal hierarchy only in (b). 



In this paper, we examine the possibility of a very long-baseline (VLBL) neutrino os- 
cillation experiment with a neutrino factory at J-PARC in resolving the neutrino mass 
hierarchy, the sign of sin 2 9 ATM — 1/2, and the degeneracy in the CP phase S MNS . The pos- 
sibility of a neutrino factory at J-PARC in Tokai is studied in Ref.j201- For the detector, 
we assume a 100 kton-level segmented water-Cerenkov calorimeter detector at L = 2, 100 

*The results in Fig ^ are obtained by using the off-axis beams JS]. Details of the J-PARC-to-HK 
analyses with the off- axis beams will be reported elsewhere 19 . 
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km away from Tokai. The distance is approximately that between Tokai and Beijing, 
where strong interests in constructing a large water- Cerenkov detector BAND (Beijing As- 
trophysics and Neutrino Detector) have been expressed [21]. The BAND detector has a 
good capability for detecting v e charged current (CC) and CC events with calorimetric 
energy measurement [2UI22]. I n this analysis, we assume that the detector is capable of 
measuring the event energy and distinguishing e from /r* 1 , but we do not require its charge 
identification capability. 

The physics prospects of VLBL oscillation experiments with a neutrino factory has been 
studied in the past by assuming that the detector can identify charges [23, and hence the 
possibility of charge identification at a large water-Cerenkov detector has been investigated 
[22] . We would like to show in this paper that even if a detector is charge blind it could 
achieve all the goals. 

In the neutrino factory, neutrinos are produced from the decay of high energy muons, 
fi + — > 77 M z/ e e + or pT — > z/ M z/ e e~. The same amount of and v e [y^ and V e ) are contained 
in the stored fi + beam. One of the most significant qualities of the neutrino factory 
is the well knowledge of neutrino fluxes. Assuming very relativistic muons, the and v e 
and V e ) fluxes from fi + (fi~) beam are expressed as 

**m = ^^jfy 1 {( 3 - 2 v) =f W - 2 v)} » ( 2a ) 
®M*e) = ^ 2 ^r2 l2 y 2 {C 1 -v)t p m(i - y)} > (2b) 

where 7 = E^/m^, y = E u /E^ with the energy of the decaying muon E^, P^ is the average 
muon polarization (P M = 1 for right-handed and P^ = — 1 for left-handed fi^), and n M is 
the number of the decaying muons. The upper (lower) sign should be taken for 
beam. 

The e-like signal from fi + beam, N e (fi + ), is given by the sum of e + from the — > V e 
appearance mode and e~ from the v e — > u e survival mode, whereas the /i-like signal, 
-^ 1 (a* + ) ) is the sum of fi + from the 17^ — > z/ M survival mode and \i~ from the v e — > 
appearance mode. The signals from the \i~ beam, N e (p,~) and N^(fi~), are obtained in 
the same way; 

Ve , (3a) 
, (3b) 
v e , (3c) 
*V (3d) 

The signals in the z-th energy bin, N{(fi + ) and NK/jT) (I = e or //), are then calculated as 
iV7(/i + ) =MN A dE v ■ P v ^ v% ■ ag c + ^ ■ P v ^ Vi ■ , (4a) 



iV e (/i H 



Ve + V 
' V„ + I 



V e + V u 



V„ 
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dE v ■ P„ ■ + $ Fe ■ ft.-p, ■ a™} , (4b) 

Ei 

respectively, where P Va ^, v and P Va ^ Vf) are the neutrino and anti-neutrino oscillation prob- 
abilities inside the earth matter, SE is the width of the energy bin, M is the mass of the 
detector, Na = 6.017 x 10 23 is the Avogadro number, and cr^ c and cr~ c are the V\ and V\ 
CC cross sections of water target [2*4*j . 

We show the dependences of the signals, Nf(p + ) and Nf(p~), on the mass hierarchy 
and sin 2 9 ATM in FigOwith 10 21 of unpolarized p + and p~~ decays, respectively, at 10 GeV. 
The left figures are the signals from the p + beam and the right ones are from the p~ beam. 
The width of the bin is taken as SE =1 GeV for E v > 2 GeV. In Fig|21 the signals for the 
normal and the inverted mass hierarchies are shown by thick and thin lines, respectively. 
We choose three values of sin 2 ATM = 0.5 (solid lines), 0.35 (dashed lines), and 0.65 (dot- 
dashed lines), where the last two values give sin 2 29 ATM = 0.91. The other parameters with 
the constant matter density p are 

Sm 2 = 3 x 10" 3 eV 2 , Sm 2 = 7 x 10~ 5 eV 2 , 

ATM ' SOL ' 

sin 2 29 SOh = 0.85 , 5 MNS = 0° , (5) 

sin 2 29 RCT = 0.06 , p = 3 g/cm 3 . 

Fig 01 shows the dependences of the expected number of events, Nf(p + ) and Nt(p~), on 
sin 2 29 RCT assuming the normal hierarchy for the same experimental setup as in FigEl We 
take three values of sin 2 26' RCT = 0.1 (dashed lines), 0.06 (solid lines), and (dotted lines). 
The other parameters are taken as eq.® and sin 2 9 ATM = 0.5. In FigOU we also show 
the number of background events by shaded bars, where In this study, we take account 
of the events coming from r pure-leptonic-decays Ni(p ± ,r — > I), neutral-current (NC) 
events, N^p^, NC), and r hadronic-decays, N[(p ± ,r — > had). The last two processes can 
contribute to e-like events, where produced 7r°'s mimic the electron shower in the detector. 

The contribution from the survival mode dominates each signal but that in N fl (p + ) and 
Nfi{^~) vanish at around E v ~ 5 GeV for sin 2 29 ATM = 1 because of the nearly maximum 
oscillation. In N e (p + ) and A^(/i~), the contributions from the appearance mode, — > V e 
and V e — > Vp, respectively, are suppressed for the normal hierarchy, because of the large 
matter effects, which suppresses Pv^Ve by a factor of 1/6 at E v ~ 5 GeV [IS], and also 
because of small anti-neutrino CC cross sections, a v Ja Ul ~ 1/2. Owing to these double 
suppressions of the anti-neutrino appearance contributions, N e (p + ) measures P Ue -^u e and 
N^(p~) measures P„ in the normal hierarchy. Because of this, N e (p + ) is very insensitive 
to the parameters sin 2 9 ATM and <5 MNS that affect the transition probabilities P v ++ Ve and 
Pn^y e . Insensitivity of N e (p + ) to sin 2 # ATM is shown clearly in FigJ21 The measurement 
of N e (p + ) hence gives us an opportunity to measure sin 2 26* RCT uniquely in the normal 
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Figure 2: Number of e- and /i-like signals for 10 21 decaying /z + (left) and yT (right) each 
at 10 GeV for the normal-hierarchy (thick lines) and the inverted hierarchy (thin lines) 
with three values of sin 2 ATM =0.5 (solid lines), 0.35 (dashed lines), and 0.65 (dot-dashed 
lines). The other parameters are taken as in eq.(JHJ). 



hierarchy. The sin 2 26 I RCT dependence of iV e (/i + ) is clearly seen in FigEl 

We find, on the other hand, that the v e and V e survival probabilities in the inverted 
hierarchy are always larger and smaller than those in the normal hierarchy, P^}°™ < 
P}?\„ and P^°™£ > P^L-n i because of the large matter effect. In fact they are almost 

£/g 'I/g i/g 'fg //g r(/g / <— ' J 

independent of sin 2 26 l RCT in the allowed region < sin 2 26> RCT <^0.1. Because P^°™ a ~ 
P l ™ and P^Z^ ~ PvLv at sin 2 26„ CT = 0, we find that the normal and inverted 
hierarchies can be distinguished by measuring iV e (/i + ) and N e (p,~), if sin 2 26 I HCT is not too 
small. The differences are quite significant for sin 2 2# RCT = 0.06, as shown in Fig|21 If the 
degeneracy in the mass hierarchy is lifted, the parameter sin 2 2# RCT can be measured with 
good accuracy by the observations of N e (fj, + ) and N e {p~), where the former is insensitive 
to sin 2 9 ATM and <5 MNS , and so is the latter for very small sin 2 29 RCT . 

The Up survival probability P v ^ Vft depends on sin 2 26* ATM in the leading term, which is 
measured uniquely by N^jjT) for the normal hierarchy owing to its small dependence on 



5 



N(H + ) 



[i beam 

e-like signal . 



23456789 10 

E v =E obs (GeV) 



NOi") 



25000 



10000 



1 


i i i i i i i i i i i i 

beam 


e-like signal 





























2 3456789 10 
E v= E obs( GeV ) 



15000 



10000- 



5000 



ji + beam 

(i-like signal 



2 3 4 5 6 7 8 9 10 

E„=E nhs (GeV) 



30000 
25000 
20000 
15000 
10000 
5000 




[i beam 
(i-like signal 



23456789 10 
E v= E obs ( GeV > 



Figure 3: sin 2 26* RCT dependences of the expected number of events. The numbers are 
calculated for eq.(0) but for sin 2 29 RCT = 0.1 (dashed lines), 0.06 (solid lines), (dotted 
lines) and sin 2 9 ATM =0.5 with the normal-hierarchy. Shaded bars show the number of 
the background events. For the r pureleptonic decay events, the horizontal scale is the 
observable energy after subtracting the final neutrino energies. 



sin 2 29 RCT as shown in FigO In Fig|2l the thick lines for the normal hierarchy show that 
the two values of sin 2 9 ATM that give sin 2 # ATM = 0.91 give almost the same prediction. This 
is because Pv R -^v is suppressed strongly by the matter effect for the normal hierarchy. On 
the other hand, in case of the inverted hierarchy, we find that the degeneracy in sin 2 29 ATM 
is resolved in N^fi^). This comes from the V e — > z7 M mode which is enhanced by the matter 
effect and receives non-negligible contribution proportional to sin 2 9 ATM ■ sin 2 9 RCT . 

The signal of iV M (/i + ) is the sum of the CC events from u e — > and — > 17 M , which 
depends on sin 2 9 ATM ■ sin 2 9 KCT and sin 2 29 ATM , respectively. Both sin 2 9 ATM and sin 2 2# ATM 
contributions can be seen in FigJ2l and sin 2 29 KCT dependence is clearly seen in FigJHl 
Because of the nearly maximum oscillation, the significant dependence on the sign of 
sin 2 # — 1/2 is seen at around E v ~ 5 GeV. 

ATM ' v 

Finally, the N e (fi~) signal in FigO shows almost linear dependence on sin 2 9 ATM for the 
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normal hierarchy, and that of Fig|3] shows a clear dependence on sin 2 2# RCT . Those are 
the results of the — > v e transition mode which is enhanced by the matter effect in the 
normal hierarchy. 

Summing up, we come up with the following strategy for resolving the degeneracies of the 
three neutrino model parameters if the neutrino mass hierarchy is normal. If sin 2 28 RCT is 
not too small, iV e (/i + ) and iV e (/i~) determine the hierarchy, and iV e (/i + ) measures sin 2 2# RCT . 
Nfj,(ii~) determines sin 2 2# ATM , and both N^(fi + ) and N e (fj,~) are sensitive to sin 2 9 ATM . We 
can foresee the sensitivity in 5 MNS as well, because the last two observables are sensitive 
to P v ^ Ve . The dependency in <5 MNS and 7r — <5 MNS may also be resolved, because their sum 
is sensitive to cos5 MNS and their difference is sensitive to sin5 MNS . A more complicated 
strategy is needed if the neutrino mass hierarchy turned out to be inverted. We will come 
back to this problem elsewhere. 

Bellow, we perform an exploratory study of the potential of such experiments by as- 
suming the following simple experimental setup: A 100 kton water-Cerenkov calorimeter 
detector at L — 2, 100 km away from a neutrino factory, which delivers 10 21 decays of 
unpolarized /i + and pT at = 10 GeV. 

We assume that the detector is capable of measuring the event energy, E Q ^ S = E^beam) — 
E u (final), with accuracy much better than 1 GeV in the interval 2 GeV < E ^ s < 10 GeV, 
and calculate the expected number of events for 8 bins with the 1 GeV width. The proba- 
bility distribution of the three neutrino parameters is then obtained from the following y 2 
function: 



x 2 = yU NKfi+r-NK^ 



' v/^ + ) obs ' { 




Sm 2 true \ 2 



l^^soi^ ] , / sin 2 2^ L -si - 2 
0.1 x Sm 2 J™ e J V °- 06 



(5m 2 fit 

_|_ I ATM 



"sol 

"sol > V 

5m 2 true \ 2 



^ L -sin 2 2g \ 2 



+ E 



a=e,fj, 



+ 



10" 4 

' e D a -l 
v 

1\ 2 



/sin 2 2fl fit - sin 2 2fl true \ 2 

_j_ I ATM ATM | 



ATM I I I ATM 

) \ 0.01 

_ > 9 



Ae 7 



<e -v 2 

C e/NC 



Ae 



:/NC 



(6) 
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Here N^/j, 1 * 1 ) ^ are calculated as 

iV:(// ± ) obs = Nid^) + NtQi* r^e)+ JV&i± r -> had) + N^, NC) , (7a) 
N;(^y hs = NiQj*) + N;(^, r - ji) , (7b) 

including contributions for the pure-leptonic r decays, and the fake contributions to iV^/i*) 
where the NC events or the v T CC events with hadronic r decays are mistaken as e ± events. 
These background contributions are calculated as 

f Ei+8E rE max 

iv;(/i + , T -+ o = mna / d£ obs / 

f do^c da cc 1 

^ M ■ *W • + ^ ■ • -j^- • B(r - ,(8a) 



Ni(fi + ,T-+ had) 




d£, 



obs 



P^ T --^ + • • "ttt 2- • B(t - had) • P e/NC , (8b) 



iV,!( / i + , NC) = MN A / 1 ' dE v • 4 C + • a* c } • P e/NC , (8c) 

and likewise for the events from decaying yT . Here E \> s is calculated by subtracting the 
energies carried away by neutrinos in r decays. The oscillation probabilities are calculated 
for a set of parameters of the three neutrino model labeled as 'true', in the normal hierarchy 
for a constant matter density of p = 3 g/cm 3 . The statistical errors of the number of events 
in each bin are simply the square root of the observed number calculated as above. 

In the calculation of the model predictions for these numbers, labeled as Ni(pi ) , we 
allow all the 6 neutrino model parameters and the matter density p to vary freely under 
the constraints given in the third and the fourth lines of eq.(|B|). For 8m 2 QL and sin 2 29 SOh , 
we set the la errors which will be achieved by the KamLAND experiment jH]. The LBL 
neutrino oscillation experiment Tokai-to-SK will determine 5r« 2 TM and sin 2 29 ATM with the 
accuracy <5(<5m 2 TM ) ^ 10~ 4 and 5(sin 2 26* ATM ) < 0.01, respectively jT3]. We assume that the 
average matter density along the baseline will be known with 0.1g/cm 3 accuracy. 

Furthermore, the following systematic effects are taken into account in our fit. Although 
we can probably neglect the flux uncertainties in the neutrino factory experiments, the 
observed numbers of events should depend on the detection efficiencies, and the expected 
numbers of events suffer from the uncertainties in the neutrino CC cross sections. We 
introduce the detection efficiency parameter and separately for 'e'-like and '//'-like 
events, respectively. For simplicity, we set = = 1 in calculating the expected number 
of events, but allow for their errors 



ef = l±Aef, e? = l±Ae?. (9) 
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Here we assume that the fiducial volume uncertainty can be accounted for as a part of 
the efficiency uncertainties. For the cross sections, we assume that the ratio a„ 
and <7p /o^ e will be determined accurately in the future and introduce the following 
uncertainty factors 



^ CC a™ a™ a cc 



- 1 ± Ae a — = - =e-=l± Ae- (10) 

( CCY fn-CCV " J-^^d / GOV rn- cc V " " A ^ 

V°V M / input V°V e J input V"77 M / input \ u V e I input 

These systematic uncertainties eq.® and eq.(fT0j) play important roles in determining the 
ultimate accuracy of the experiments when the statistical errors get smaller, and appear 
in the fifth line of eq.(jOJ). Below, we show our results for Ae® = Ae® = Ae G v = Ae^ = 
and 2 %, so that the impacts of reducing these systematic uncertainties can be inferred. 

We also take account of the systematic errors for the uncertainty in the branching 
fractions of r leptonic-decays as 

B(r^eu e u T ) _ B(j l±£ ^^ 



B(T -> eV e V T ) input B(T -> llffiPr) input 

with B(t — > ez7 e z/ T ) input = 0.178 and B[r — > /i77 M i/ T ) inpu t = 0.174 [26J, and Ae T = 0.1. 
For the NC events, we set the mean value of the e/NC misidentification probability as 
(-^e/Nc)input = 0.25% by using the estimations from the K2K experiment 0. The uncer- 
tainty factors are then introduced as 

N< " r e/Nc = 1 ± Ae e/NC , (12) 



(Pe/Nc) 



input 



with Ae e /Nc = 0.1. Contributions from these systematics are accounted for as the last two 
terms of eq.fjHJ). 

We show in FigEJx 2 ^ as a function of sin 2 29^ when the event numbers calculated for 
the normal hierarchy are analyzed by assuming the inverted hierarchy. The results shown 
in Fig0] (a) are obtained when the detecting efficiencies and the cross sections are known 
exactly, As? = Ae? = AeZ = Ae£ = 0. Because of P™°™ < PJ*T U and P^°™± > Pi n ^ , 
the normal and inverted hierarchies can be distinguished by measuring N e (fi + ) and N e (/jT). 
Normal hierarchy can be established at 5a level when sin 2 26^™° J> 0.01, and at 3a level when 
sin 2 2$^ ;> 0.004. However, the differences of the e-like signals between the normal and 
the inverted hierarchy is reduced to ~ 1% when sin 2 2# RCT = 0.004. It is hence necessary 
to measure the cross section of u e and V e CC events with an accuracy better than 1%. In 
Fig El (b), we show the results when we account for 2 % uncertainties in the cross sections 
and the efficiencies, Ae? = Ae? = Ael = ^ e v = 0.02. The hierarchy discrimination 
power of the experiment is reduced significantly, when the 5(3) a discrimination is now 
possible only when sin 2 20*™ £0.1(0.008). 
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Figure 4: Xmm as a function of the fitting parameter of sin 29^ by assuming the inverted 
hierarchy. The resluts are for 10 21 fi + and fi~ decays for a 100 kton water target at 
L = 2, 100 km. The input data are calculated for the 'true' values in eq.(JT3J) by assuming 
the normal hierarchy The 12 sets of the input data are labeled by sin 2 2f9* r "° = 0.01 , 0.004 
and for (a), 0.02 ,0.008 and for (b), with 5^ g = 0° (solid lines), 90° (dotted lines), 180° 
(dashed lines), 270° (dot-dashed lines). Aef = Ae£ = Ae a v = Ae£ = for (a), and 
Aef = Ae° = Ae° v = Ae£ = 0.02 for (b). 



The input data for Fig0] (a) and FigQ] (b) are generated by assuming the normal hier- 
archy with the following 'true' values : 



sin 2 9 true 

ATM 
,2 o/Qtrue 
SOL 



sin z 29 



0.5, 
0.85 



5m 2 true = 

ATM 
,2 true 



5m„ 



SOL 



3 x 10" d eV z , 
= 7 x 10~ 5 eV 2 



sin 2 20!™° = 0.01, 0.004, (for (a)), 0.02, 0.008, (for (b)) 



rtrue 

MNS 



0°, 90°, 180°, 270° 



(13a) 
(13b) 
(13c) 
(13d) 



The 12 sets of the input data are labeled as sin 2 29*™° = 0.01 , 0.004 and in Figg] (a) and 
0.02,0.008 and in FigH(b) with 5^™ = 0° (solid lines), 90° (dotted lines), 180° (dashed 
lines), 270° (dot-dashed lines). We then fit the data by assuming the inverted hierarchy by 
allowing all the parameters to vary freely, and obtain the x 2 ■ values plotted in the figures. 

We show in Fig|S] Xmin as a function of sin 2 9^ for (a) Aef = Ae^ = Ae a v = Ae^ = 0, 
and (b) Aef = AeJ = Ae a v = Ae£ = 0.02. Input data are calculated for sin 2 6*^ = 
0.35, 0.4, and 0.45 (sin 2 20^™ = 0.91,0.96, and 0.99, respectively) with three values of 
sin 2 20^ = 0.06 (solid lines), 0.02 [0.08] (dotted lines) and 0.01 [0.1] (dashed lines) for 
(a) [(b)] and = 0° with the normal hierarchy. The values of the other parameters are 
taken as in eq. (|T5j) . The Xmm function is found by varying the fitting parameters within 
the normal hierarchy. We see that each Xmm h as two dips at the value of sin 2 which 
give the same sin 2 26* ATM . The results in Fig|5] (a) show that we can resolve the degeneracy 
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Figure 5: xLin as a function of the fitting parameters sin 2 9^ TM for (a) As® = As® = 
Ae a v = Ae£ = and (b) Aef = Ae^ = Ae a v = Ae£ = 0.02. The input data of the events 
are calculated for sin 2 9 tIue =0.35, 0.4, and 0.45 with sin 2 2# true = 0.06 (solid lines), 0.02 

ATM ' ' RCT v ' ' 

[0.08] (dottted lines) and 0.01 [0.1] (dashed lines) for (a) [(b)] and = 0°. The other 
input values are the same as in eq.(JT3J). 



at 3a level when sin 2 ^ = 0.4 (0.45) for sin 2 2^ ^0.02 (0.04). The degeneracy in 
sin 2 ATM can be resolved by the u e — > mode in A^(/i + ) and — > u e mode in N e (fi~) 
whose leading terms are proportional to sin 2 9 ATM sin 2 9 BCT , when sin 2 29 BCT and sin 2 20 ATM 
are determined by the survival mode in N e (fi + ) and iV^(/i~), respectively. In Fig|5] (b) 
where we take account of 2 % systematic errors in the CC cross sections and the detection 
efficiencies, the degeneracy in the sign of sin 2 ATM — 1/2 can be resolved at 3a level when 
in 2 0true = o.4 f or s i n 2 2# true ^0.06. In the case of sin 2 9 tIue = 0.45, the determination of 

ATM RCT ~ ATM ' 

the sign is possible at 3a level only when sin 2 26^™ ~ 0.1. 

In FigEl we show allowed regions in the plane of sin 2 29^ and <5*j* when sin 2 26^™° = 
0.02 , 0.06 , and 0.1 with = 0° (upper left), 90° (upper right), 180° (lower left), and 270° 
(lower right) in the normal hierarchy and the other parameters are eq. ()13a|) and eq. (j!3b|) . 
The errors in the efficiencies and the cross section are set as Ae® = Ae® = Ae a v = Ae^ = 
0.02. In each figure, the input parameter points are shown by solid-circles. The normal 
hierarchy is assumed in the fitting. The regions where Xmin < 1, 4, and 9 are depicted by 
solid, dashed, and dotted boundaries, respectively. 

The figures for = 90° and 8%^ = 270° show that the CP phase 5 MNS can be 
constrained locally around the 'true' points. We can discriminate the maximal CP violation 
cases 5^ e s = 90° , 270° from the CP conserving cases <5 MNS = 0° , 180° at 3a level for 
sin 2 26 , J^ ^0.02. The sensitivity is based essentially on the T violating difference between 
Pve-^v^ and P„ -+ Ve . Because the matter effects enhance both of these transitions whereas 
suppressed strongly the transitions 77 M — > V e and V e — > z/ M , A^(/x + ) is sensitive to P V£ ^ Vli 



11 



360 



270 















• 






\ "J: 


























































































































































































































































& 






/ d 








• ''/, 


.YV. 




) 


12 


0.04 0. 


)6 0.08 





1 0. 



=180 



90 



0.02 0.04 0.06 0.08 

sin228 RCT 



0.1 





























































































































• 










» 






1 























































































































































0.12 



360 



270 



90 



0.02 0.04 0.06 0.08 0.1 0.12 

sin220 RCT 



Figure 6: Allowed regions in the plane of sin 2 26 I RCT and 6 MNS for sin 2 20*^ = 0.02 , 0.06 , 
and 0.1 with = 0° (upper left), 90° (upper right), 180° (lower left), and 270° (lower 
right). The errors in the efficiencies and the cross section are taken to Aef = Ae^ = 
Ae° = Ae^ = 0.02. In each figure, the input parameter points are shown by solid- circles. 
The regions where Xmin < 1, 4, and 9 are depicted by solid, dashed, and dotted boundaries, 
respectively. 



and N e (n ) is sensitive to P v _> Ve . 

For <5 truc = 0° and 5 trne = 180°, we can constrain 5..,„ within ±45° accuracy for 

MNS MNS MNb J 

sin 2 2# truc = 0.06 and 0.1, while the fake solutions appear for sin 2 2# tme = 0.02 at 3a 

RCT ' RCT 

level. This is essentially because both v e — > and — > u e transitions have terms 
proportional to sin 2 29 BCT cos 5 MNS ; N e (fi + ) measures sin 2 29 KCT uniquely, and the sum of 
N^(fi + ) and N e (/j,~) is sensitive to sin 2 26 I RCT cos 5 MNS while the difference is sensitive to 
sin 2 29 RCT sin <5 MNS . In our previous study of J- PARC superbeam and Hyper-Kamiokande 
[T5] , it was difficult to resolve the degeneracy which gives same values for sin 2 26' RCT cos <5 MNS . 

In this letter, we examine the capability of the VLBL neutrino experiment with a neu- 
trino factory at J-PARC in Tokai Village to resolve degeneracies in the neutrino oscillation 
parameters. We find that a large segmented water- Cerenkov calorimeter detector placed 
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at a few thousand km away, such as a proposed BAND detector [21] in Beijing, can be very 
effective in solving the problem. The charge identification capability at the detector is not 
required. This is mainly because the e ± events from the fj, + beam measure the v e — > v e 
survival rate rather accurately while the e ± events from the \x~ beam are more sensitive 
to the Up — > v e transition rate, for the normal hierarchy, where the earth matter effects 
strongly suppress the — > V e and V e — > transitions. On the other hand, if the hierarchy 
is inverted, the — > v e and v e — > transitions are suppressed and the e ± events from 
fi + beam are sensitive to the F M — > V e transitions and those from fi~ beam are sensitive 
to V e — > V e . Because of the strong matter effects in the <-> v e and <-> V e transitions 
processes, we can resolve not only the neutrino mass hierarchy but also the degeneracies 
in the sign of sin 2 ATM — 1/2 and those between 5 MNS and 180° — S MNS . 

The following numerical results are found at 3 a level for 10 21 decaying unpolarized 
/i + and /i~ at 10 GeV when we take into account 2% uncertainties in the detection ef- 
ficiencies and the v e and CC cross section measurements and assume that the energy 
resolution is significantly smaller than 1 GeV. The mass hierarchy can be determined for 
sin 2 2# RCT ^0.008, which is remarkably smaller than the corresponding limit in the VLBL 
experiment with the conventional beams from J-PARC [T3"j . When sin 2 2# RCT ^0.06, the 
degeneracy in the sign of sin 2 6* — 1/2 can be lifted for sin 2 29 = 0.96. The CP 

& J & ATM ' ATM 

violating phase 6 MNS can be uniquely constrained for sin 2 20 J>0.02 if its true value is 
around 90° or 270°, while if it is around 0° or 180°, the mirror solution at 180° — S MNS can 
be excluded for sin 2 26' RCT ^0.03. It is essential to have both /i + and \i~ beams to obtain 
the above results, but it is not necessary for the detector to have the charge identification 
capability. 
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